Electric field control and optical signature of entanglement in quantum dot molecules 
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The degree of entanglement of an electron with a hole in a vertically coupled self-assembled dot 
molecule is shown to be tunable by an external electric field. Using atomistic pseudopotential calculations 
followed by a configuration interaction many-body treatment of correlations, we calculate the electronic 
states, degree of entanglement and optical absorption. We offer a novel way to spectroscopically detect 
the magnitude of electric field needed to maximize the entanglement. 

PACS numbers: 78.67.Hc, 73.21. La, 03.67.Mn 



c3 



-a 
c 
o 

(N 
> 

oo 

(N 
O 
in 
o 

-I— > 
c3 



i 

C 

O 

o 



X 



A pair of quantum dots or "quantum dot molecule" 
(QDM) occupied by two electrons lup] or by an electron- 
hole pair j3l have been offered dUUas a basis for quan- 
tum computing. The fundamental requirement for such 
quantum algorithm is the availability of entangled states 
and the ability to entangle and disentangle the quantum bits 
(qubits). In the context of a dot molecule, an entangled 
electron-hole pair can be represented by the the maximally 
entangled Bell state e T h T + e^/i^, where e and h stand 
for the electron and the hole (the two qubits) and T and 
B for their localizations in top or bottom dot. The origi- 
nal proposal 131 and subsequent experiments 0, Eh for 
entangled electron-hole pairs in QDMs promised a high 
degree of entanglement [3] based on analysis via simple 
models. However, later theoretical work showed that 
electron-hole entanglement is generally low in such cases 
and develops a sharp maximum only at a specific interdot 
separation that critically depends on the size difference of 
the two dots. Unfortunately, it has proven to be difficult 
to experimentally control so precisely the interdot distance 
and the size difference of the two dots. The question we 
address here is whether the degree of entanglement can be 
maximized by other means, more accessible experimen- 
tally than a variation of the interdot separation. We pro- 
pose and quantify theoretically that it is possible to tune 
and control the degree of entanglement by applyingan ex- 
ternal electric field in the growth direction II7L IsLM [Toll . 
The use of electric field has been demonstrated in quantum 
dots 17ll8l.l9l. il OB and very rec ently in a single quantum dot 
molecules by Krenner et al. ill ill . We predict that, while 
the entanglement at zero field is generally low (35% for 
our case), it can reach a high value (75% in our case) at 
a specific electric field Fs ma x (-5.4 kV/cm in our case). 
Moreover, precisely at this field the first two exciton lines 
merge, giving a well defined spectroscopic signature of the 
point of maximum entanglement. 

In order to obtain reliable results for the correlated ex- 
citon states, it is of foremost importance to accurately 
account for the multi-band character of the hole states 
and for the correct strain dependence in the coupling re- 
gion (between the dots). We have thus solved the perti- 
nent Schrodinger equation atomistically, in a multi-band 
fashion. We use the Hamiltonian H = — 1/2V 2 + 



Yla n v a( r ~ Rn) + Kso + \e\Fz under an external elec- 
tric field F applied in [001] (z) direction. The atomistic 
pseudopotentials v a of atom of type a and the non-local 
spin-orbit potentia l Va n are fit to reproduce InAs and GaAs 
bulk properties ll6l ll2ll . The atomic positions {H n } are ob- 
tained by minimi zing the atomistic strain energy (via va- 
lence force field 113ll ) for a given shape and size of the 
dots. Our quantum dots have a truncated cone shape (12 
nm base and 2 nm height) with a composition ranging from 
pure InAs at the top to Ino.5Gao.5As at the base, as de- 
termined in Ref. The single-particle Hamiltonian is 
diagonalized in a basis ^ = ^ k A nik (f) nik of pseudopo- 
tential Bloch functions <j) n ik as outlined in Ref. [l4j, thus 
permitting coupling of various Bloch states. Correlations 
are treated via a many-body expansion in Slater determi- 
nants Jl5l fl6ll where the electrons not included dynami- 
cally are represented by a model screening of the Coulomb 
and (long and short range) exchange [17]. The entangle- 
ment is calculated according to the von Neumann entropy 
of entanglement |lj| UjJ]. 

The bonding (6) and antibonding (a) electron molecular 
levels of a dot molecule will be denoted as E a , E b . For 
an idealized (mostly unrealistic) symmetric case the lowest 
energy molecular orbitals (MOs) develop from single-dot 
electron states and e B located on the bottom (B) and 
top (T) dots: 

ip[E b ] = -L(e T + e B ); i/j[E a ] = -J=(e T - e B ), (1) 

and similarly for the holes H a , H b . As shown previ- 
ously (a |2(J], in reality, because of strain and random- 
alloy fluctuations, one does not have a symmetric bonding- 
antibonding behavior even if the dot molecule is made 
of identical (but non-spherical) dots. This is seen in 
Fig-Hwhere both electron and hole molecular orbital wave 
functions are shown for zero electric field F = 0. We 
see that the (lighter- mass) electrons tunnel between dots, 
forming bonding-antibonding states as in Eq. Q, but the 
(heavier-mass) holes remain localized on the top (bottom) 
dot for the bonding (antibonding) MO H b (H a ). The 
single -particle molecular orbital energy levels are shown 
in Fig. |3 As we apply an electric field the molecular lev- 



2 





E = kV/cm 


E= 3.6 kV/cm 


E=9 kV/cm 




t — u KV/cm 


t — 3.o KV/CIT1 


t — y KV/cm 































LU 








I 
























































LLJ 








X 

























FIG. 1: Square of the first two electron and first two hole wave 
functions of an InGaAs/GaAs quantum dot as a function of elec- 
tric field at an interdot separation d=l .9 nm. 




-15 -10 -5 5 10 15 

Electric Field in [001] direction (kV/cm) 



FIG. 2: Single particle electron (upper panel) and hole (lower 
panel) eigenvalues as a function of electric field for d=l .9 nm 
with respect to the GaAs valence band maximum. We indicate 
the localization of the MOs on top and bottom dot by ex, e B , 

h T , h B . 



els that develop from the single-dot orbitals exhibit anti- 
crossing. We have indicated in Fig.|2]fhe major character 
of the molecular states E a , E b , H a , H b in terms of the lo- 
calization on individual dots (er, e B , hx and h B ) using 
the calculated MO wave functions of Fig. \l\ We see that 
for holes at positive fields ip(H a ) ~ h B and ip(H b ) — Ha 
while for electrons ip(E a ) ~ and ip(E b ) ~ e B . The 
opposite is true for negative fields. Thus, by applying an 
electric field we can tune the localization of the MO's and, 
for instance, compensate for size, composition or shape dif- 
ferences of both dots. We will see that this tuning of local- 
ization will also control the degree of entanglement. 

There are four transitions between the four molecular 
levels shown as vertical arrows in Fig. |2] Their single- 
particle transition energies e 1 ^ (differences between the en- 
ergies from Fig.|2J are given in Fig. [3 a) and show maxima 
and minima vs field. We note in Fig. |3ja) the character 
of the four transitions in terms of localization on single- 
dot orbitals. We see that at high fields, the lowest- and 
highest-energy transitions involve different dots: for exam- 
ple E b H b is eshr at positive fields, and erh B at negative 
fields. Thus, the corresponding dipole transitions are ex- 
pected to be weak ("dark states"). In contrast the second 
and third transitions at high fields involve the same dots: 



(a) Single-Particle Transition Energies 
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(d) Correlated Excitonic Energies 
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FIG. 3: (Color) Transition energies for d=l .9 nm as a function 
of electric field in different approximations: (a) Single particle 
transition energies e g (b) Direct electron-hole Coulomb matrix 
elements between MOs J e h . (c) Transition energies including 
electron-hole Coulomb interaction: e 1 ^ + J**!, but without cor- 
relation effects, (d) Final correlated exciton results. The bidot 
products (eTfiT, etc..) are given whenever the MO states (given 
on the right) are strongly dominated by one of these products. 



for example E b H a is e B h B at positive fields and e T h T at 
negative field. Thus, the corresponding dipole transitions 
are expected to be large ("bright states"). 

The single -particle approximation underlying Fig. |3a) 
is valid only in the case of large fields, where e-h Coulomb 
effects are small compared to the field-driven variation in 
the single-particle levels. At these large fields, the excitons 
are pure determinants with localization on either e T h T or 
e T h B or e B h T or e B h B and therefore show no entangle- 
ment. We will next see that in the interesting region of 
electric fields, e — h Coulomb interactions are crucial. In 
Figure|3b) we show the calculated electron-hole Coulomb 
interaction 



Je,h [i - j] 



e(r a ,r b )\r a - r b \ 



dr a dr b 
(2) 



3 



between the MOs (i = E„ or E b ) and (j = H a or H b ) 
using the model of Resta lfl7ll for the screening e. This 
interaction J e h (Fig. 0b)) shows the reverse behavior vs 
field compared with the MO energies e g vs field [Fig.[3ja)]. 
For example, whereas J e ^ h [E b H b ] is maximal (less neg- 
ative) at large positive or negative fields, and minimal at 
intermediate fields, the MO band gap e g [E b H b ] is minimal 
at large positive or negative fields and maximal at inter- 
mediate fields. Not surprisingly, when one calculates the 



Coulomb-corrected excitonic transition energy e g 



J e,h 



between the molecular states i and j [Fig. 0c)] one sees 
a partial cancellation for the two low-energy transitions, 
E b H b and E b H a , leading to a weak dependence of the tran- 
sition energy on field. In contrast, inspection of Figs. 0a) 
and 0b) for the two highest-energy transitions shows that 
the field-dependence of e*'- 7 ' and J*'^ reinforce each other, 

so the excitonic gap e 1 ^ + J l e '\ [Fig. 0c)] has an amplified 
dependence on field. We conclude that the combination of 
e g and J e h brings the lowest-energy transitions closer to 
each other, while pushing the two higher-energy transitions 
apart. This will affect the correlation coupling between the 
MO's, as seen next. 

The Coulomb-corrected excitonic transition energies 
£*'-' + J\ \ neglect the interactions between the different 
configurations, i.e., the states from Fig. 0c) are not allow 
to interact. This interaction is included in the next step via a 
configuration interaction rcpfilli^i calculation in which 
we include all Coulomb and exchange integrals from the 
first four electron and first four hole states (including spin). 
The results are shown in Fig. 0d) as a function of elec- 
tric field. The lowest energy transitions (excitons |1) and 
1 2)) have a very weak dependence on field, similarly to the 
case without correlations [Fig. 0c)]. In contrast to the per- 
turbative approach of Fig. 0c) the states do not cross but 
anticrossing at -5.4 kV/cm, as expected from interacting 
states. Similarly, |3) and |4) anticross at+3.6kV/cmbutin 
a more abrupt fashion. 

To understand the correlated CI results we next ana- 
lyze |1), 1 2), 1 3), 1 4) by decomposing the correlated exci- 
tonic states into sums of products of single-dot states e^/iy, 
enhr, exhs and e^/i^ called "bidot products" (see Ref. 
Oil for details). The results of this decomposition are given 
in Fig. 0a)(b) for states |1) and |2) as a function of the 
electric field. Fig. 0c) shows the degree of entanglement 
calculated for the correlated CI wave functions using the 
von Neumann formula 1(4 11811 S = — Tr p e log 2 p e , where 
p e is the reduced density matrices of the electron. When a 
state is made solely of a single bidot product such as eshr 
it is unentangled, but when it is made of a coherent su- 
perposition, such as esliB ± erhr, it might be entangled. 
In the case of very strong positive fields (larger than 20 
kV/cm), state 1 1) is purely e#^T (unentangled), as the field 
pulls the electron to the bottom dot and the hole to the top 
dot. The excitons |2) and |3) are already, at moderate pos- 
itive fields, purely eshs an d &rhr, respectively, and re- 
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FIG. 4: Occupation probabilities given as the bidot products 
e^hx or exliB or e^/ir and eghg as a function of the electric 
field for the first two exciton states |1) (a) and j2) (b) for d=7.9 
nm. (c) Entropy of entanglement of the first two excitons. 



main unentangled states at large fields. The entanglement 
of |1) and 1 2) [Fig. 0c)] reaches its maximum of around 
75% at F Smax = -5.4 kV/cm. At this field the exciton states 
|1) and 1 2) are mainly composed of erhr ± ^b^b con- 
figurations, as shown in Fig.0a)(b) [23]. The calculated 
excitonic states of Fig. 0d) are now used to calculate the 
absorption spectra in Fig.0 In Fig.0the oscillator strength 
is plotted for different values of the electric field as a func- 
tion of the transition energy. The plot shows a total of four 
transitions marked with |1), |2), |3), |4). Since the tran- 
sitions 1 3) and |4) are weak, their position is marked by 
solid and open arrows respectively. The transitions 1 1) and 
|2) (|3) and |4)) have a weak (strong) dependence on field 
and show an anticrossings at E Sms , x (Eg PT ). The spectra 
show that at £ Sma) . the lowest energy exciton 1 1) becomes 
dark and progressively gains oscillator strength away from 
the anti-crossing. The point of merging of |1) and |2) at 
the field E Smax reflects a "resonant conditions" with maxi- 
mum entanglement (Fig.0. The distinct spectroscopic sig- 
nal of the anticrossing (where the lowest line progressively 
looses oscillator strength) occurs at the point of maximum 
entanglement and, we suggest, can give experimentalists a 
simple way to control a delicate quantity such as entangle- 
ment. 
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FIG. 5: Oscillator strength of the first four transitions (first 16 
transitions including spin) as a function of electric field for d=7.9 
nm. 
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FIG. 6: Formation of the highly entangled excitonic states at the 
critical field F=-5.4 kV/cm. (a) Single particle electron and hole 
levels in the dot-localized basis, (b) Simple differences between 
the single particle electron and hole energies from (a), (c) Adding 
electron-hole direct Coulomb interaction to (b). (d) Adding elec- 
tron and hole hopping the the levels from (c) 



The physics underlying the "resonant condition", that 
produces the high degree of entanglement, is revealed in 
Fig. using the more intuitive basis of dot-localized or- 
bitals (as opposed to MO's). To obtain dot-localized or- 
bitals, we rotate the MO's (that are delocalized over both 
dots) until the on-site Coulomb interaction is maximized 
l2lll . This procedure yields single -particle energies of the 
dot-localized orbitals, denoted as e T , e B , hx and h B in 
Fig- Eta). Fig. 0a) shows that at Fs max the energies of the 
dot-localized electron and hole orbitals are separated by 10 



and 7 meV respectively with a gap of 1.3 eV. Generally, 
the energy separation between e T and e b an d between /i T 
and h B can reflect size, composition or shape differences 
of the two dots. These differences can be tuned by the elec- 
tric field. Fig. 0b) shows the energies of simple products 
of these electron and hole states. They show two closely 
spaced levels (3 meV apart) in the center of the spectrum 
and two states 7 meV lower and higher in energy. These 
energies are different from the MO's energies of Fig.QJa) 
that are combinations of exhx, e^s, &b^t an d eg/tg 
at this field. In the next step, in Fig. 0c), the two-body 
Coulomb attraction is taken into account and lowers the 
e T h T and states in such a way that they are about 

14 meV below the e^s state. This is the consequence of 
a weak e — h binding for the dissociated excitons exhs 
and eshr- Notably, the simple products e T h T and eshs 
are energetically nearly degenerate at this level, this is the 
"resonant condition" mentioned above. In the last step, in 
Fig.0d), the excitons |1) and |2) are now created by in- 
cluding the effects of electron and hole hopping that effec- 
tively produce correlated states. The excitons |1) and |2) 
are now a bonding- and antibonding-like combination of 
the energetically degenerate e^hx and eshs states. The 
excitons |1) and |2) are now split by a small energy of 
less than 1 meV. This small splitting is conceptually very 
similar to the Davydov splitting lE^l in molecular crystals. 
The analysis also reveals that 1 1) is anti-symmetric (e T \iT - 
e^s) an d therefore optically dark while |2) is symmetric 
(e^/iT + eshs) and optically bright. The high symme- 
try of these states (purely symmetric and anti-symmetric) 
leads to the high degree of entanglement. Any deviations 
from Fsmax will l ea d to a less symmetric combinations as 
-^(ae T h T + (3e B h B ) with a ^ f3 with lower entangle- 
ment and smaller oscillator strength. 

In conclusion, we showed that the degree of electron- 
hole entanglement in coupled quantum dots can be tuned 
by an external electric field and that the point of maximum 
entanglement can be identified by measuring the photolu- 
minescence spectra, observing the merging of two peaks. 
This opens new ways for experimentalists to identify the 
electric field needed to achieve maximum entanglement in 
specific dot molecules. We finally analyzed the nature of 
the excitons and revealed the interplay of single particle 
effects, direct Coulomb binding and electron and hole hop- 
ping on the many body levels. We described how these 
effects conspire to yield a highly entangled state. 

This work was supported by the US Department of En- 
ergy (SC/BES) under the LAB 03-17 initiative. 



[1] D. P. DiVincenzo, Science 270, 255 (1995). 
[2] D. Loss and D. P. DiVincenzo, Phys. Rev. A 57, 120 (1998). 
[3] M. Bayer, P. Hawrylak, K. Hinzer, S. Fafard, M. Korkusin- 
ski, Z. R. Wasilewski, O. Stern, and A. Forchell, Science 



5 



291, 451 (2001). 
[4] K. Hinzer, M. Bayer, J. P. McCaffrey, P. Hawrylak, M. Ko- 

rkusinski, O. Stern, Z. R. Wasilewski, S. Fafard, and 

A. Forchel, Phys. Stat. Sol. (b) 224, 385 (2001). 
[5] M. Bayer, G. Ortner, A. Larionov, V. Timofeev, A. Forchel, 

P. Hawrylak, K. Hinzer, M. Korkusinski, S. Fafard, and 

Z. Wasilewski, Physica E 12, 900 (2002). 
[6] G. Bester, J. Shumway, and A. Zunger, Phys. Rev. Lett. 93, 

47401 (2004). 

[7] P. W. Fry, I. E. Itskevich, D. J. Mowbray, M. S. Skolnick, 
J. J. Finley, J. A. Barker, E. P. O'Reilly, L. R. Wilson, I. A. 
Larkin, P. A. Maksym, et al., Phys. Rev. Lett. 84, 733 
(2000). 

[8] I. Shtrichman, C. Metzner, B. D. Gerardot, W. Y. Schoen- 

feld, and P. M. Petroff, Phys. Rev. B 65 (2002). 
[9] M. Sugisaki, H. W. Ren, S. V. Nair, K. Nishi, and Y. Ma- 

sumoto, Physical Review B 66 (2002). 
[10] B. Alen, F. Bickel, K. Karrai, R. Warburton, and P. Petroff, 

Appl. Phys. Lett. 83, 2235 (2003). 
[11] H.J. Krenner, M. Sabathil, E.C. Clark, A. Kress, D. Schuh, 

M. Bichler, G. Abstreiter, and J.J. Finley, Phys. Rev. Lett. 

94, 057402 (2005). 
[12] A. J. Williamson, L.-W. Wang, and A. Zunger, Phys. Rev. B 

62, 12963 (2000). 
[13] P. N. Keating, Phys. Rev 145, 637 (1966). 



[14] L.-W. Wang and A. Zunger, Phys. Rev. B 59, 15806 (1999). 
[15] A. Szabo and N. S. Ostlund, Modern Quantum Chemistry 

(McGraw-Hill, New York, 1989). 
[16] A. Franceschetti, H. Fu, L.-W. Wang, and A. Zunger, Phys. 

Rev. B 60, 1819 (1999). 
[17] R. Resta, Phys. Rev. B 16, 2717 (1977). 
[18] C. Bennett, H. Bernstein, S. Popescu, and B. Shumacher, 

Phys. Rev. A 53, 2046 (1996). 
[19] G. Bester, A. Zunger, and J. Shumway, Phys. Rev. B 71, 

075325 (2005). 

[20] L. He, G. Bester, and A. Zunger, cond-mat/ 0503492, 
(2005). 

[21] C. Edmiston and K. Ruedenberg, Rev. Mod. Phys. 35, 457 
(1963). 

[22] D. P. Craig and S. H. Walmsley, Excitons in molecular crys- 
tals (W. A. Benjamin, INC., New York, 1968). 

[23] In order for this e-h entanglement to be useful in a quan- 
tum computation scheme, the qubits need to be physically 
separated to be addressed individually. One possibility for 
such a separation seems to be the use of in plane electric 
fields that would coherently drive the electron and hole to 
different neighboring quantum dots. This would constitute 
the next experimental challenge 



